Bilayer graphene heterostructure with hBN: 
Interplay between misalignment, interlayer 
cn asymmetry, and trigonal warping 
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We show that a slightly incommensurate hexagonal substrate (e.g., hexagonal boron ni- 
tride) applied to one side of bilayer graphene (BLG) with Bernal stacking inflicts inversion 
^ symmetry breaking and opens gaps in the spectrum of minibands created by the periodic moire 



pattern. This contrasts the behaviour of moire minibands in monolayer graphene, where well- 



q resolved minibands are connected by Dirac points. 
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The heterostructures of graphene with other hexagonal layered crystals or crystals with hexag- 

•i-H 

onal symmetry facets feature moire patterns which are the result of incommensurability of the pe- 
riods of the two two-dimensional lattices, or their misalignment. Since the period of the moire pat- 
tern is longer for a pair of crystals with a closer size of lattice constants and better aligned principal 
crystallographic axes, the long period moire superlattices are characteristic for graphene/hexagonal 
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boron nitride (hBN) heterostructures with a small misalignment angle between the two honey- 
comb lattices. Such heterostructures have recently been created by transferring graphene onto 
hBN. 1_4 The influence of hexagonal moire patterns on Dirac electrons in monolayer graphene 
(MLG) has been studied in detail, 1_8 both using specific microscopic models and phenomenolog- 
ically. Three possible types of moire miniband structures on the conduction/valence band sides of 
graphene's spectrum have emerged from the theories: 2 ~ 8 sometimes, spectra without a distinct sep- 
aration between the lowest and other minibands; quite exceptionally, the first miniband separated 
from the next band by a triplet of secondary Dirac points (sDPs) in each of the graphene valleys K 
and K'; more generically, a single sDP at the edge of the first miniband in each valley. Also, the 
signatures of the miniband formation have been observed experimentally in the tunneling density 
of states 2 and magnetotransport characteristics 3 ' 4 in MLG/hBN heterostructures. 

In this Letter, we analyze the characteristic moire miniband features in heterostructures of bi- 
layer graphene (BLG) with highly oriented and almost commensurate, hexagonal crystals, such as 
hBN, recently created and investigated using magneto-transport measurements by Dean et al. 9 We 
find that, in contrast to monolayers, the electronic spectrum of BLG on hBN is most likely to ex- 
hibit global gaps between the first moire miniband and the rest of the spectrum (on either valence 
or conduction band side, but sometimes in both bands), or have the bands strongly overlapping 
with each other, but almost never has Dirac points at the miniband edge. This behaviour is pre- 
scribed by the substrate creating a moire perturbation only for one layer of BLG, thus breaking the 
inversion symmetry of the moire superlattice. This point becomes apparent from the parametric 
space diagrams in Fig. 1, where the regions of the parameter space with gapped spectra are painted 
in red with an overlapping (non-resolved) bands - left transparent. Differences between the two 
diagrams corresponding to different misalignment angles 6 arise from the interplay between the 
orientation of the supercell Brillouin zone (sBZ) and interlayer hopping in BLG. 

The above conclusions are derived using a phenomenological approach 8 which involves the 
description of the long-range moire superlattice using a Dirac-type model for graphene electrons, 
where we include the dominant harmonics in the moire perturbation and perform an exhaustive 
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Figure 1: Parameter space (uo,ui,uj) used to classify characteristic behaviour of moire miniband 
in highly oriented BLG on almost commensurate substrate: in red, we paint regions where BLG 
spectrum has a gap separating the first miniband in the valence band from the rest of the spectrum, 
regions with overlapping (unresolved) bands are transparent and blue/green mark the degenerate 
conditions for the appearance of isolated secondary Dirac points at the first miniband edge, like 
in monolayer graphene. 2 ' 6 " 8 A similar parametric plot describing the minibands on the conduction 
band side can be obtained by inversion in the U0-U3 plane. 

numerical analysis in order to characterise the miniband behaviour over a broad range of the pa- 
rameter space. For BLG placed on top of a substrate with hexagonal symmetry and the lattice 
constant cis = (l + 8)a larger by 8 than that of graphene, (a = 2.46A, and for the case of hBN, 
5=1.8% ), the lattice mismatch, together with a possible misalignment of the two lattices given 
by the angle 0, lead to a periodic structure which can be described using a set of reciprocal lattice 
vectors, 
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where Ry stands for anticlockwise rotation by angle (p and b = \b n \~A=-\/S 2 + 2 . 4 Note that this 
set both rotates by 0(0) and changes its size as a function of 0, Fig. 2. The electrons in BLG/hBN 
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Figure 2: Hexagonal Brillouin zone (BZ) of BLG with the valley coordinate system p x ,p y and the 
electronic bands in the vicinity of the K valley, together with the zoom in on the supercell Brillouin 
zone (sBZ) where we mark its own symmetry points K and /i. 
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written in the basis 10 of the Bloch states on sublattices {<j> (Ai), (B2) , (A2), (#1 ) } in the K val- 
ley and {(j)(B2),—(j>(Ai),(j)(B\),—(j)(A2)} in K' , where indices 1/2 mark the bottom/top layers and 
the substrate directly acts on the electrons in the bottom layer. We also use fi = 1 and employ two 
sets of Pauli matrices a,, o = (o x , o y ), and %u acting in the sublattice and valley space, respectively. 
The first matrix on the left-hand side of Eq. (??) follows from the tight-binding description of 
BLG, ! 1 with v— 10 6 m/s resulting from the interlayer hopping in graphene. The 'vertical' interlayer 
hopping between the closest neighbours is taken into account by parameter 71 =0.38eV, 12 and 
the skew next- neighbour interlayer hopping, 73, is taken into account by the 'trigonal warping 



term', with V3 ~ O.lv. This results in the four-band BLG model, with two bands degenerate 

at the points K and K' ', and two bands split away by ±71. 

The next two terms on the right-hand side of Eq. (??), scaled using energy scale vb and 



parametrised using dimensionless wo,i,2,3> describe the effect of the substrate on BLG through 
its coupling with carbon orbitals in the closest, bottom layer. The parameter uq characterizes the 
magnitude of the perturbing potential, uj, captures the local asymmetry between the Ai and B\ 
sublattices, whereas u\ and u^ introduce modulation of the in-plane hops for the electrons travel- 
ling within the bottom layer. As opposed to the monolayer case, 8 the U2 term cannot be completely 
gauged away but its effect on the band structure is generated only via the interplay with the trigonal 
warping term, and hence is small, so that we neglect it in further discussions. Note that, since we 
scale all energies using the scale vb set by the periodicity of the moire pattern, the size of dimen- 
sionless parameters m, in Eq. (??) would be larger for smaller angle 6 for the same pair of BLG 
and substrate crystal. 

For the sake of a systematic comparison with the low-energy 'two-band' model of free-standing 
BLG, 11 applicable at the energy scale £<Cyi, we use a Schrieffer- Wolff transformation 14 and 
project the four-band Hamiltonian onto the low-energy bands, reducing it to an effective two-band 
Hamiltonian, 
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The applicability of the simplified Hamiltonian to the description of, at least, the first moire mini- 
band in the BLG spectrum requires that j\ > 2vb. For a perfectly aligned BLG/hBN heterosystem, 
we estimate that \ ~ ^^-—jw^w. = 1-08, which suggests that a quantitative description of moire 
minibands in BLG/hBN requires the use of the four-band Hamiltonian, Eq. (??). 

From this point, the analysis of moire minibands in BLG consists in a systematic zone folding 
and numerical diagonalization of the four-band Hamiltonian H. Due to the time-inversion of the 
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(e) 0=0.597°, 0=-3O°, vb = 1.087j 
w = -0.06, Mi = -0.06, w 3 = 0.03 
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Figure 3: (a-e) Moire miniband spectra (drawn within the rhomboidal sBZ) and density of states 
(DoS) for 6 =0 and 6 =0.597° and several points in the (wo, wi, W3) parameter space as marked in 
Fig. 1 . Also, an example of the spectrum for large misalignment angle is shown in (f). In cases 
(b,d,f), a sizeable gap is visible between the conduction and valence gap of BLG, opened due to 
the interlayer asymmetry generated by the moire perturbation of the substrate applied only to one 
layer. 
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15 



H( P )=a y Ty[H(-p)] T GyTy, 



electronic spectra in the two valleys are related, £k+ p = £> , and in the following, we only discuss 
electronic spectra in valley K. Moreover, using the commutation properties of the Pauli matrices, 



(a) 6=0, 0=0, vb = 0.93y 1 

M = 0.2, u 1 = u 3 = 0,5 = -0.041eV 




0.5 




I ' 25 


^^ 


ffi o 








QJ 




S-0.25 


^. 


^—^ 


-0.5 





DoS [arb. units] 
(b) 0=0.597°, ^=-30°, i# = 1.08^ 

uo = -0.06, «i = -0.06, u 3 = 0.03, 5 = -0.056eV 

0.5F 




DoS [arb. units] 



Figure 4: Examples of the miniband spectra and DoS for structures doped up to four holes per 
moire supercell using a single side gate. By inspection, one can find that constant in space in- 
terlayer asymmetry, leading to a sizeable gap between valence and conduction bands of BLG, 
does not affect the occurrence of a gap or the overlap between the first and next minibands in the 
superlattice spectrum. 

one finds that for the band structure within the same valley, 
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This relation partly simplifies the exploration of the parameter space (wo,«i,«3) anu we use it to 
relate the parametric maps for the minibands properties in the valence and conduction bands in 
graphene. 

Figure 3 illustrates characteristic behaviour of the moire miniband on the valence band side in 
BLG/hBN for various choices of perturbation parameters and misalignment angles: (a) overlap- 
ping minibands characteristic for the transparent part of the parameter space in Fig. 1; (b) gapped 
spectrum at the edge of the first moire miniband characteristic for the red-painted part of the pa- 
rameter space; (c) rare occurrence of a sDP, at the boundary between the above two characteristic 
regimes. This behavior is shown for both perfectly aligned graphene and hBN lattices (0 = 0) in 



Fig. 3(a-c), and for 6 = 0.597° in Fig. 3(d,e). Increasing the misalignment angle leads both to a 
shorter period of the moire pattern (reflected by smaller size of dimensionless parameters w,) and 
to a rotation of the sBZ, to = -3O° for = 0.597° and = -6O° for = 1.836° (vfc= 1.89ft). 
We find that this rotation of the moire pattern interplays with the anisotropy of the warping term 
v^o-p in Eq. (??), additionally affecting the occurrence of gapful and gapless miniband spectra. 
Note that for all angles, there is a gap at zero energy, 

A ps 6vbuo(iiQ + u 3 ), (3) 

between the conduction and valence bands. Opening the gap at the neutrality point of BLG, in the 
third order of perturbation theory, is analogous to opening a small gap by the inversion symmetry 
breaking part of the moire perturbation in monolayer graphene, but with the difference that the 
moire pattern itself does not have to be inversion asymmetric: the inversion symmetry is broken 
by the application of substrate perturbation to only one layer out of two. For BLG sandwiched 
between two hBN layers, the asymmetry will be still present, if the misalignment angles are differ- 
ent, and both the gap A and the low-energy dispersion will be determined by that hBN layer which 
has a smaller misalignment angle. This is because, for the same crystalline substrate, the size of 
parameters u\ is smaller for larger vb, and the result of Eq. (??) suggests scaling A<* d~ 2 . 

A clear separation between minibands in the spectrum of a periodically modulated two-dimensional 
semiconductor results in a change of sign of the Hall coefficient as a function of doping 3 ' 9 as the 
band filling passes the critical values of four electrons per moire supercell and the chemical poten- 
tial jumps from the lower miniband to higher (here we take into account spin/valley degeneracy). 
Such a behaviour has been observed in magneto-transport measurements by Dean et al., 9 in a 
similar fashion to monolayer graphene/hBN heterostructures. 3 Since this happens at a relatively 
high carrier density, 3 ' 9 for which the required large gate voltage introduces an additional interlayer 
asymmetry, a gap A between the conduction and valence band is further enhanced, which can af- 
fect spectroscopic characteristics of the system. This purely electrostatic effect. 11 ' 16 is specific 



to devices with a single-side gate, whereas the spectrum in devices with compensating top and 
bottom gates would not be affected by doping. To take into account this gate-induced interlayer 
asymmetry, we add an additional term to Eq. (??), 
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and estimate A for the carrier density of four electrons per moire supercell using the mean field 
theory, 16 ' 17 and find that its presence has almost no effect on the type of the miniband spectrum at 
the edge of the first miniband (gapped or overlapping minibands), as one can conclude by compar- 
ing spectra in Fig. 4 with those in Fig. 3. Therefore, the conclusion on the occurence of the two 
types of miniband spectra over various parts of the parameter space (uq,u\,u^) in BLG remains 
virtually the same as illustrated in Fig. 1. Since the systematic study of the moire miniband spectra 
in BLG shows that the parameter space is dominated by the occurrence of two types of spectra: 
with a gap separating first and the next miniband, or, alternatively, strongly overlapping minibands 
(as in weakly modulated semiconductor structures in GaAs/AlGaAS 18 " 21 ), here, we argue that a 
change in sign of the Hall coefficient in BLG/hBN should be expected to accompany an activation 
behaviour of conductivity o xx . This is in contrast to monolayer graphene, where the hBN sub- 
strate induces secondary Dirac points at the first miniband edge and thus leads to a finite value 
in the conductivity maximum at the carier density corresponding to four holes per moire super- 
cell (see Ref. (9)). Surprisingly, the latter feature of BLG (activation behavior of conductivity at 
the superlattice-induced resistivity maximum) has not been clearly identified in the experiment in 
Ref. (9). 
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